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. Abstract. In this paper we use a design theoretical approach to construct new, previously 

unknown complex Hadamard matrices. Our methods generalize and extend the earlier 
results of [TJ, [TU] and offer a theoretical explanation for the existence of some sporadic 
f T | ' examples of complex Hadamard matrices in the existing literature. As it is increasingly 

difficult to distinguish inequivalent matrices from each other, we propose a new invariant, 
the fingerprint of complex Hadamard matrices. As a side result, we refute a conjecture of 
Koukouvinos et al. on (n — 8) x (n — 8) minors of real Hadamard matrices [15] . 
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1. Introduction 

A complex Hadamard matrix if is a square n x n matrix with arbitrary unimodular 
■ entries with complex orthogonal rows, i.e. HH* = nl. They are the natural generalization 
of real Hadamard matrices. Constructions of complex Hadamard matrices are motivated by 
their various applications in quantum information theory [18], harmonic analysis [8j, |16j . 
operator theory [12], and combinatorics [20]. Let us recall that Hadamard matrices H and 
K are equivalent if H = P1D1KD2P2 holds for some permutational matrices P\,Pi and 
unitary diagonal matrices D\,D2. An Hadamard matrix with its first row and first column 
consisting only of Is is said to be normalized. 

The pragmatic example for complex Hadamard matrices are the (rescaled) Fourier ma- 
trices F n , and more generally the Butson type Hadamard matrices, all of whose entries are 
some fixed m th roots of unity [3J. As the concept itself is still a discrete generalization of 
real Hadamard matrices, many aspects of the real theory can be successfully applied to them 
obtaining interesting structural-, existence- and non-existence type theorems (9], [19] . An 
infinite family of "exotic" (i.e. non Butson type) matrices, the circulant complex Hadamard 
matrices, were discovered in the early 90s by the pioneering work of Bjorck [1] and the 
follow-up papers [7] and [TO] . Recently it was pointed out by Dhja that the so-called gener- 
alized Kronecker product construction leads to matrices with free parameters, and therefore 
examples of non Butson type complex Hadamard matrices in composite dimensions [5J- The 
only examples of parametric families of prime orders are Petrescu's biunitaries |llj . 

The outline of the paper is as follows. In Sections 2-4 we use design theoretical methods 
to construct new examples of complex Hadamard matrices. Besides rediscovering some well- 
known matrices of small orders, we give a new and systematic proof for the existence of 
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non-standard circulant complex Hadamard matrices of prime orders [7J, [10], and we present 
new, previously unknown families of complex Hadamard matrices as well. In particular, 
we show that the matrices ?7 15 , Vi 5 and Wqa are new to the literature. As it is increasingly 
difficult to distinguish inequivalent complex Hadamard matrices from each other, in Section 5 
we introduce a new invariant, the fingerprint of complex Hadamard matrices. By combining 
Craigen's paper jl] with the theory developed in Section 5 we show, as a side result, that 
the (n — 8) x (n — 8) minors of real Hadamard matrices of order n cannot take the values 
k ■ 2 7 ■ n n / 2 ~ & for k £ {28, 29, 30, 31} disproving a conjecture of Koukouvinos et al. [13]. 

Through the paper we shall use the standard notations for well-known matrices such as 
F n , C n , P n etc. The description of these matrices is available in the online version of the 
Tadej-Zyczkowski catalogue of complex Hadamard matrices [15], [17] . Notations C/4 m _i, 
V4 m _i, W^ m+ i)A and W^ m+ \)B are to be introduced in this paper. 

2. Block designs and complex Hadamard matrices 

In this section we consider symmetric 2-designs with the usual parameters v, k, A. We 
always identify a block design B with its incidence matrix B. In this way we simply consider 
0-1 matrices of order v, with k Is in every row and column satisfying BB* = (k — X)I + A J, 
where J is the all 1 matrix. To exclude trivial cases, we suppose, as usual, that A < k < v . 

It is well-known that the algebraic object behind real Hadamard matrices is the 2- (4m — 
1, 2m — 1, m — 1) Hadamard design. The combinatorial nature and the high level of internal 
symmetry of the design ensures that it can be trivially transformed to a real Hadamard 
matrix of order 4m, by simply exchanging every with —1, and padding the obtained matrix 
with a full row and column of Is. It is natural to ask whether we can construct complex 
Hadamard matrices in a similar fashion as well. In the following we give a characterization 
of complex Hadamard matrices, composed of two different entries. The main concept of this 
section relies on the following 

Definition 2.1. We say that a block design B induces a complex Hadamard matrix if after 
exchanging every in B with a fixed complex number a of modulus one, the obtained matrix 
is complex Hadamard. We say that this is an induced complex Hadamard matrix. 

Recall that a complex Hadamard matrix is regular, if the absolute value of the sum of 
the entries in each row is constant. The following two lemmata indicate that there is a one- 
to-one correspondence between block designs satisfying ([T|) and regular complex Hadamard 
matrices, composed of two different entries. 

Lemma 2.2. A 2-(v, k, A) design induces a complex Hadamard matrix if and only if 

(1) v-Vv <k< v + Vv^ 
K J 2 - - 2 

Proof. Consider a block design B. After exchanging every with a, the orthogonality relation 
between any two rows of it reads 

(2) X + 2(k- \)ft[a]+v + \-2k = 0. 

From this, by using A = kik — l)/(f — 1), one can express 3? [a] in terms of v , k, and as \a\ = 1 
should hold, we have 

< 3 > «w-i-$^-i. 
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After rearranging the desired inequality follows. 

Clearly, induced Hadamard matrices are regular. We have the following 



□ 



Lemma 2.3. Any regular complex Hadamard matrix H of order v, composed of entries 
{l,a}, corresponds to a 2-(v,k,X) block design. 

Proof. If a = —1, then H is real, and the corresponding object is the Menon design. Oth- 
erwise, if a 7^ —1, then an easy (but tedious) analysis shows that the regularity condition 
implies that every rows and columns contains the same number of Is, say k. Now consider 
any two rows of H, and suppose they share a common 1 in exactly A coordinates. The same 
orthogonality relation (T5]) implies that for these two rows we have A = k + v/(23t[a] — 2). 
Therefore the value of A is independent of the choice of the rows. □ 

Observe that the parameters of the Hadamard design satisfy the inequality (JT|). We have 
the following 

Theorem 2.4. Suppose that we have a 2-(4m — 1, 2m — 1, m — 1) Hadamard design, repre- 
sented by an incidence matrix U. Then, after replacing every with 



(4) 



-1 



1 . V4m - 1 



2m 

in U we obtain a complex Hadamard matrix U, 



± i 



2m 



4m— 1 • 



It is well-known that the Paley-I construction gives rise to circulant Hadamard designs 
when p = 3 (4) is a prime. Using this fact the authors of [10] constructed circulant com- 
plex Hadamard matrices of prime orders, without exploring the possibilities of constructing 
Hadamard matrices of composite orders as well. In particular, we have the following 

Corollary 2.5 (Munemasa-Watatani, [TO])- For every prime p = 3 (4),p > 7 there exists a 
circulant complex Hadamard matrix of order p, inequivalent to the Fourier matrix F p . 

Remark 2.6. It is an important open problem in the theory of operator algebras to decide 
whether there exist infinitely many inequivalent complex Hadamard matrices of prime orders 



Applying Theorem 12.41 for m = 1,2,3 we get the matrices F 3} C 7 a,Cib and Cua,Cub 
respectively, while for m = 4 we obtain from the 5 inequivalent Hadamard designs at least 5 
new, previously unknown complex Hadamard matrices of order 15. We exhibit a particular 
example obtained from the five fold tensor product of F 2 : 



(5) 



7 . 
-- +i 
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So far we have shown that regular complex Hadamard matrices, composed of two different 
entries, correspond to block designs. It is natural to ask whether there are examples of non- 
regular complex Hadamard matrices, composed of two different entries as well. Somewhat 
surprisingly, such matrices exist in the real case only. 

Lemma 2.7. Suppose that we have a complex Hadamard matrix H, composed of two different 
entries, say {l,a}. If a ^ —1, then H is regular. 

Proof. Consider the inner product of any two rows of H: 

(6) A + Ba + Ca + D = 0, 

where A,B,C,D are integral numbers describing in how many coordinates meet the pairs 
(1,1), (I, a), (a, 1) and (a, a), respectively, i.e. A denotes the number of columns sharing 
a common 1 in the rows considered, etc. As a is non-real, we have B = C, and therefore 
A + B = A + C . This shows that the number of Is is the same in every two rows of H, hence 
H is regular. □ 

Remark 2.8. Let n := k — A be the order of a symmetric design. It is well known, that 
An — 1 < v. On the other hand, equation ([2]) implies v < An. Therefore induced complex 
Hadamard matrices correspond to the Hadamard- and Menon designs onlyQ 



3. Conference matrices redux 

In the previous section we have described the rediscovery of the so called cyclic p-roots of 
"index 2" type matrices of orders p = 3 (4). It is quite natural then to try to look at the 
underlying structure of the circulant complex Hadamard matrices when p = 1 (4) as well. To 
do this, we invoke an other design theoretical object, the conference matrices. Recall, that 
a conference matrix C is a square matrix with 0s on the main diagonal and ±1 otherwise 
satisfying CC* — (n — 1)1. A conference matrix is normalized, if all the nonzero entries in 
the first row and column are 1. We are interested in symmetric conference matrices of orders 
4m + 2. The main result of this section is the following 

Theorem 3.1. Given any normalized, symmetric conference matrix C of order Am + 2 > 6 
one can construct a complex Hadamard matrix in the following way: discard the first row 
and column of C replace the 0s with 1, replace the off-diagonal Is with c and replace the 
off- diagonal —Is with c where c is an unimodular complex number with 



(7) RH 1±^™±I. 

Am Am 

This procedure give rise to complex Hadamard matrices Wu m +UA an d Wu m +i)B) depending 
on the sign of ([7]). 

Proof. Consider any normalized, symmetric conference matrix of order Am + 2, and after 
neglecting its first row and column replace the diagonal 0s to 1, replace the off-diagonal Is 
and —Is to unimodular indeterminates x and y respectively, obtaining a matrix W. Our 
aim is to show that by setting x = c, y = c we get a complex Hadamard matrix. To do this, 
consider any two rows of W, and by pre- and post multiplying it by the same permutational 
matrix we can suppose that the considered rows are the first two. Now we have two essentially 
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different cases, depending on either Wip = x or W12 = V hold. In the first case the scalar 
product of the rows reads 

(8) m{x\ + 2m^[xy] + 2m - 1 = 0, 
while the second case leads to the equation 

(9) 23%] + 2mU[xy] + 2m - 1 = 0. 

Solving the system of equations (jSJ)-© yields the desired result. □ 



Remark 3.2. The reader might amuse himself by checking that starting from m 
W / (4m+i)A and W(4 m+ i)B are inequivalent. 



2, matrices 



The construction corresponding to circulant conference matrices was discovered earlier in 
[7]. The authors of that paper, however, did not explore the possibilities of constructing 
complex Hadamard matrices in composite dimensions. We have the following 

Corollary 3.3 (de la Harpe- Jones, [7]). For every prime p = 1 (4),p > 13 there exists a 
circulant complex Hadamard matrix of order p, inequivalent to the Fourier matrix F p . 

The existence of matrices F 5 , Ci^b follows from Theorem 12.41 However, already for 
m = 2 we obtain a new complex Hadamard matrix of order 9. With choosing the positive 
sign in (jTJ), we have 



(10) 



9A 



1 

4 



15 



It is standard to show that this matrix is not included in the Tadej-Zyczkowski catalogue. 
In particular, it is inequivalent from Ng. The other choice of the sign in (JTJ would lead to a 
Butson type Hadamard matrix, composed of third roots of unity. 



4. A GENERALIZATION 

In Section [2] starting from a single orthogonality condition, namely equation §2§ we con- 
structed complex Hadamard matrices with two different entries. In Section [3] we started 
from a less "regular" combinatorial object, and by introducing 3 different entries we ob- 
tained Hadamard matrices with two essentially different orthogonality relations. In this 
section we start from block designs again, but we sacrifice some of the internal symmetries 
of the design in order to get still a feasible number of orthogonality equations. We have the 
following 

Theorem 4.1. Given any normalized, symmetric, real Hadamard matrix H of order km > 8, 
one can construct a complex Hadamard matrix in the following way: discard the first row 
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and column of H , replace all off-diagonal — Is with b and replace all diagonal Is with 
where b is the following unimodular complex number: 



-b, 



-l 



± i 



. \jAm — 5 



2m - 2 2m - 2 

This procedure gives rise to complex Hadamard matrices V^ m ^\. 



Proof. The proof is similar to the calculations carried out during the proof of Theorem 13.11 
but slightly longer, as in this case we have 6 essentially different orthogonality equations. 
After replacing the diagonal Is with x, the diagonal —Is with y and finally, the off-diagonal 
— Is with z, the arising orthogonality equations to be satisfied (up to conjugation) are the 
following: 

2m^[z] + 2&[x] + 2m - 3 = 0, (m > 3), 



(12) 
(13) 
(14) 
(15) 
(16) 
(17) 



2(m - + 2$t[xz] + 2m - 1 = 0, 

2mM[z] + x + y-z + 2m-2 = 0, 
2m^t[z] + x~z + zy — ~z + 2m — 2 = 0, 
2(m - i)$t[z) + 29%] + 2m - 1 = 0, 

= 0. 



2mM[z] + 29ft [yz] + 2m - 

Solving the system of equations f|T2|) - f|T7|) is straightforward. One can obtain either {x, y, z} = 
{1, a, a} corresponding to the construction described in Theorem 12 .41 or {x, y, z} = {— 6, —1, b}, 
as desired. □ 

For m = 2, 3 the construction described in Theorem 14.11 leads to the rediscovery of Pe- 
trescu's matrix P7 (see Example 15.4ft and one of Nicoara's sporadic matrices Nua- Also, for 
m = 4 we can construct at least one new, previously unknown complex Hadamard matrix 
of order 15, as follows 
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Again, it is straightforward to see that this matrix is not included in the Tadej-Zyczkowski 
catalogue. These small order examples shows that Theorem 12.41 and Theorem 14.11 lead to 
inequivalent complex Hadamard matrices. We shall emphasize this fact in the following 
section. 

Remark 4.2. One might try to investigate skew-symmetric matrices as well in a similar 
fashion as described in Theorem 14. 1L However, in that case nothing else than Theorem 12.41 
can be obtained. 
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5. Equivalence of complex Hadamard matrices 

From time to time new complex Hadamard matrices appear in the literature, and it is 
increasingly difficult to realize that a "newly discovered" matrix is indeed inequivalent to all 
existing ones. The methods available to deal with the real case cannot be applied directly 
to complex matrices, as they either based on the combinatorial nature of real Hadamard 
matrices (and therefore cannot be generalized at all), or become computationally expensive 
after proper modifications. In this respect we seek for a new invariant which might help, in at 
least the small order cases, to quickly identify classes of inequivalent Hadamard matrices from 
each other. The basic method for detecting inequivalence of complex Hadamard matrices 
was the idea of Haagerup's [6], who considered 2x2 submatrices and a complex valued 
function defined on them. More precisely, the following set A associated with a complex 
Hadamard matrix [H]ij = hij of order n 

(19) A(H) = {hijhkihuhkj : i,j,k,l = l,2,...,n} 

is invariant under equivalence. As an illustration, we show the following 

Proposition 5.1. Complex Hadamard matrices £/4 m _i andV± m _\ are inequivalent ] e or m > 2. 

Proof. Suppose that 4m — 1 is not a square. Normalize the matrix U^m-i and observe 
that there should be a matrix element u G Q[iy/4m — 1] \ Q[i]. Clearly u G A(f/ 4m _i) but 
u ^ A(V4 m _i). Otherwise, if 4m — 1 is a square, then 4m — 5 is not, and we can repeat, 
mutatis mutandis, the same proof starting with V^m-i- D 

On the one hand this invariant conveniently distinguish Butson Hadamard matrices com- 
posed from different roots of unity from each other, on the other hand, however, it cannot 
detect inequivalence between matrices composed from the same roots of unity. In particular, 
A(H) = {±1} for every real Hadamard matrix (of order at least 2). The features of complex 
Hadamard matrices characterized by the Haagerup set are "local" in the sense that only 
2x2 submatrices are considered. In order to capture at least some of the global properties 
of a complex Hadamard matrix, we shall consider higher order submatrices as well. The 
downside of our approach is that we do not introduce any fancy complex function defined on 
the submatrices, but we use the absolute value of the determinant instead. Let us introduce 
the following invariant of ours. 

Definition 5.2. Let H be a complex Hadamard matrix of order n > 4. For fix 2 < d < \n/2\ 
let 1(d) be an index set and let us denote by Vt(d) and rrii(d),i G 1(d) the absolute values of 
the d x d minors of H , and their multiplicity, respectively. Then, the following ordered set 

(20) $(H) := {{( Vi (d), mi (d)) : i G 1(d)} : d = 2, . . . , [n/2\ } 

is the fingerprint associated to H. Note that the pair (vi(d) , rrii(d)) is ordered as well. 

Naturally, for every d and for every i,j G 1(d), i ^ j it follows that v,i(d) ^ Vj(d). Note 
that 1(d) simply "counts" how many different values are taken by the absolute values of the 
d x d minors of H. 

As permuting and rephasing the rows and columns of a matrix does not affect the absolute 
values of its minors, it is straightforward to realize that $ is a class function. 

Proposition 5.3. The fingerprint of a complex Hadamard matrix is invariant under the 
usual equivalence. 
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Example 5.4 (Petrescu's matrix, [llj). As it was explained earlier Theorem \4-l\ implies the 
existence of the following matrix 
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For d = 2 it is easily seen that the absolute values of the 2x2 minors take four different 
values only. Therefore the cardinality of 1(2) is 4, and we are free to set 1(2) = {1, 2, 3, 4}. 
These four values with their respective multiplicities read v± = 0,mi = 54; v 2 = 1,^2 = 
114; t> 3 = 2,m 3 = 96 and t> 4 = a/3, 1714 = 177. For d = 3 the absolute values of the 3x3 
minors take 12 different values already. Hence, we can set the index set 1(3) = {1,2,..., 12}. 
The fingerprint, associated to Pj reads 

(22) $(P 7 ) = { {(0, 54), (1, 114), (2, 96), (V3, 177)} , {(0, 60), (1, 36), (2, 108), (3, 210), 

(4, 110), (V3, 162), (2\/3, 216), (3^3, 14), (V7, 111), (v^, 54), (Vl~9, 36), (V21, 108)}} . 

With the notations used in Definition 15.21 it is clear that the sets corresponding to the 
cases d = 1 and d = n are simply {(l,n 2 )} and {(n n//2 , 1)} respectively, and therefore they 
are omitted from the definition of $. Somewhat less obvious, however, why the remaining 
cases \n/2\ < d < n are excluded as well. This is explained in the following 

Proposition 5.5. For 1 < d < n — 1 with the notations of Definition \5.2\ we have 

(23) {(vi(n -d),m i (n-d)) :i e I(n - d)} = {(n n/2 -\(d) , m^d)) : i e 1(d)}. 

The proof of Proposition 15.51 easily follows from the following two linear algebraic lem- 
mata. One suspects after realizing (|23|) that in Hadamard matrices there is a one-to-one 
correspondence between the minors of size d and n — d, which is indeed the case. The first 
lemma is the well-known generalized matrix determinant 

Lemma 5.6. Suppose that A is an invertible n x n matrix, U, V are n x m matrices. Then 

(24) det(A + UV*) = det(J + V*A~ 1 U)det(A). 

Despite our best efforts, we were unable to find any references to the following 
Lemma 5.7. Given any unitary matrix 



(25) 



U 



A B 
C D 



with blocks A, B, C, D where A and D are not necessarily of the same size square matrices, 
then we have |det(A)| = |det(D)|. 
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Proof. As U is unitary, we have 

(26) AA* + BB* = I 

(27) CC* + DD* = I 

(28) AC* + BD* = O 

(29) A*B + C*D = O, 

where the last equation follows from the fact that U* is unitary as well. If both A and D 
are singular, then we are done. Otherwise we can suppose that, for example, A is invertible. 
Hence from (1291 we have B = — (A*)~ 1 C*D, and then, by plugging into (1251) we get 

(30) A*AC*C = C*DD*C. 

Then, starting from equation ( J26l) we use repeatedly Lemma 15.61 to obtain 

(31) |det(A)| 2 = det(J - BB*) = det(J - B*B) = det(J - D*C(A*A)~ 1 C*D) = 

= det{A*A - C*DD*C)/det{A*A), 
which, by formula (13"Uj) equals to 

(32) det{A*A - A*AC*C)/det{A*A) = det(I - C*C) = det(J - CC*) = \det{D)\ 2 . 

□ 

Proposition 15.51 has some striking applications in the determination of minors of real 
Hadamard matrices. For example, one might ask whether a 6 x 6 ±1 matrix with maximum 
determinant of 5 • 2 5 can be embedded as a submatrix into a 8 x 8 real Hadamard matrix 
This is easily seen to be impossible, as by formula (1231 Vi(6) = 8 8 / 2 2 v { (2),i G 1(2) 
should hold. As v j(2) G {0, 2} we conclude that fj(6) can assume the values or 4 • 2 5 only. 
In particular, it is enough to study the distribution of minors up to size n/2 in Hadamard 
matrices. Apparently the authors of [13] were unaware of this fact, and although they 
observed that there is some connection between j x j and (n — j) x (n — j) minors (for small 
j), they conjectured that the n — 8 minors can take the values k ■ 2 7 ■ n n ^ 2 ~ 8 , k = 1, 2, . . . , 32. 
Again, by Proposition 15.51 we have Vi(n — 8) = n n ^ 2 ~ 8 v j(8), i £ 1(8), however, a result of 
Craigen [1] shows that Vi(8)/2 7 <£ {28,29,30,31}. In particular, k G {28,29,30,31} in the 
conjecture above is impossible. 

We conclude our paper with the following remark: the classification of all cyclic proofs 
of index 3, and the related circulant complex Hadamard matrices have been completed very 
recently [2J. It would be interesting to see what kind of combinatorial object lies behind 
those matrices, if any. 
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